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The non-equilibrium response of a superconductor to pair-breaking radiation 
measured over a broad frequency band 
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We have measured the absorption of terahertz radiation in a BCS superconductor over a broad range of 
frequencies from 200 GHz to 1.1 THz, using a broadband antenna-lens system and a tantalum microwave 
resonator. From low frequencies, the response of the resonator rises rapidly to a maximum at the gap edge 
of the superconductor. From there on the response drops to half the maximum response at twice the pair¬ 
breaking energy. At higher frequencies, the response rises again due to trapping of pair-breaking phonons in 
the superconductor. In practice this is the first measurement of the frequency dependence of the quasiparticle 
creation efficiency due to pair-breaking in a superconductor. The efficiency, calculated from the different non¬ 
equilibrium quasiparticle distribution functions at each frequency, is in agreement with the measurements. 


In a superconductor at low temperature, most of the 
electrons are bound in Cooper pairs. These pairs can 
be broken into quasiparticles by absorbing photons with 
an energy larger than the binding energy. This mecha¬ 
nism is frequently used to detect submillimetre and tera¬ 
hertz radiation using conventional superconductors such 
as aluminium. Pair-breaking detectors are usually as¬ 
sumed to measure the number of quasiparticles created 
by the absorbed radiation. The observable that measures 
the number of quasiparticles varies from the complex 
conductivity for microwave kinetic inductance detectors 1 
(MKIDs), the current through a tunnel junction 2 to the 
capacitance of a small superconducting island 3 . These 
observables are mainly sensitive to quasiparticles with 
an energy close to the gap energy of the superconductor, 
A. The working principle of these detectors is usually ex¬ 
plained in terms of an effective number of quasiparticles 
which is maintained by a balance between the radiation 
power and electron-phonon interaction (recombination) 4 . 
To convert the power (P) into a number of quasiparti¬ 
cles (N qp ), the quasiparticle creation efficiency r\ p b is in¬ 
troduced, which compares the actual N qp with the maxi¬ 
mum possible N qp when all created quasiparticles would 
have an energy A. Since Cooper pairs have a binding 
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energy of 2A, a photon with an energy in between 2A 
and 4A can still only create two quasiparticles. The rest 
of the energy is lost through electron-phonon scattering, 
hence r] p b < 1. For higher energies rj p b depends on the 
phonon trapping factor, which determines whether high 
energy phonons are directly lost or can break an addi¬ 
tional pair, rjpb is therefore not an efficiency in the sense 
that photons are lost, but it reduces the detector respon- 
sivity. MKIDs are superconducting microwave resonators 
which sense the number of quasiparticles through the 
complex conductivity of the superconductor. The phase 
response (9) of such a resonator can be approximated by 

9 oc oc 5N qp oc ri pb P, (1) 

02 

where 02 is the imaginary part of the complex conductiv¬ 
ity. For the last proportionality we assume to be in the 
linear regime where the quasiparticle recombination life¬ 
time does not change significantly upon a change in N qp . 
N q p is dominated by background power and 5N qp <C N qp . 
T] p b , and hence the detector response, is dependent on the 
frequency of the absorbed photons, even at constant ab¬ 
sorbed power. 

On a microscopic level, the pair-breaking radiation 
leads to injection of quasiparticles at very specific 
energies 5 . Together with electron-phonon interaction 
(scattering and recombination) 6 a non-equilibrium, non- 
thermal quasiparticle energy distribution f(E) is formed, 
which determines the response to pair-breaking radiation 
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FIG. 1. (a) Schematic of the detector. The detector chip is 
fabricated on a SiN membrane and glued to an elliptical lens, 
leaving a small vacuum gap between the antenna and the lens- 
dielectric. (b) The design of the detector chip. The antenna 
slot, coupled to a microwave resonator, is fabricated on a SiN 
membrane (light grey square). The resonator is capacitively 
coupled to a microwave readout line. 


as recently shown in Ref. 7. For microwave resonators 
this is reflected in the explicit dependence of 02 on f(E) 8 : 
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where cfn is the normal state conductivity, Ti Planck’s 
constant and uj the microwave frequency. r] p b is thus an 
attempt to capture all information contained in f(E) in 
a single number, to allow for an effective quasiparticle 
number approach as given by Eq. 1. 

Here we present the first measurement of 77 ^ over a 
broad range in frequencies close to the superconducting 
gap (350 - 1100 GHz). A Ta MKID is used as the detector 
in a Fourier transform spectrometer (FTS) to measure 
the frequency dependence of the response. The measured 
response curve of the detector can be well explained by 
a frequency dependent 77 ^, caused by a different non¬ 
equilibrium f(E) calculated for different pair-breaking 
frequencies. 

From an applied point of view MKIDs 9 are considered 
promising detectors for large arrays due to the intrinsic 
ease of multiplexing their readout. MKIDs are photon 
noise limited for various frequencies 10 15 . The level of 
experimental detail that has now been achieved 13,16,17 
calls for a more detailed understanding of the absorption 
of radiation. An important gap in this understanding is 
a measurement of 77 ^. r] p b determines key parameters: 
the responsivity of the detector, the recombination noise 
level in the photon-noise limited regime 10 and the sen¬ 
sitivity in the generation-recombination noise dominated 
limit 18 . The common number used for 77 p b is 0.57 for 
all signal frequencies, which was derived for the tempo¬ 
ral relaxation of very high energy excitations which first 
create a photo-electron 19,20 , an approach which is not 
applicable for frequencies close to the gap. 

Previous studies of the absorption of radiation in su¬ 


perconductors have either measured f(E) directly with 
tunnel-junctions 21,22 , but only with a single-frequency 
optical laser, or measured the absorption over a broad 
band with a bolometer 23,24 , which is insensitive to the 
non-equilibrium effects that determine r] p b- To measure 
ripb over a broad frequency band, a known and relatively 
constant radiation power over a broad frequency band 
is required. Secondly, we need the absorption of all of 
that power at all frequencies within the volume of the 
detector to exclude the effect of frequency dependent 
absorption 23 . We therefore use a particularly wideband 
lens-antenna system, which is based on the leaky-wave 
antenna 25 and shown in Fig. 1. It consists of a 30 fi m 
wide, 4 mm long slot, etched in a 200 nm thick Ta film 
with a resistivity of 6.7 /if2cm, which is sputter deposited 
onto a 3 /im thick SiN membrane using a 6 nm Nb seed 
layer. A spacer chip, placed in between the Ta and the Si 
lens ensures a 35 /im vacuum gap between the metal layer 
and the Si lens, which is crucial to get a high directiv¬ 
ity of the antenna over a broad frequency band 25 . This 
lens-antenna was demonstrated to have very clean beam- 
patterns over the frequency range 300-900 GHz. The 
antenna launches the signal as a travelling wave into a 
coplanar waveguide (CPW) with a central strip of 3.5 
/im and slots of 3 /im wide, which length is designed to 
make a quarter wavelength resonator at 4.6571 GHz (the 
MKID detector). An extensive discussion of the design, 
fabrication process and beampattern measurements can 
be found in Ref. 26. 

The detector is cooled down in a 3 He/ 4 He cryostat to 
a bath temperature of 320 mK. The cryostat has optical 
access through a window, Goretex infrared blockers at 77 
K and 4 K, a 1.1 THz lowpass filter at 4 K. The Michel- 
son Fourier transform spectrometer (FTS) consists of a 
globar source at 2000 °C, a fixed and a movable mirror 
and a mylar beamsplitter. To eliminate absorption lines 
due to water, the FTS is placed in vacuum. The MKID 
itself is the detector in this setup. The phase response 
of the detector was measured as a function of the mirror 
distance. The phase response is linear in power, which is 
verified using the response to a full rotation of a polariser 
in a separate measurement (i.e. the last proportionality 
in Eq. 1 is valid). The Fourier transform of the inter- 
ferogram, corrected for the frequency dependence of the 
filters and beamsplitter (see Supplementary Figure SI 27 ), 
is shown in Fig. 2 as black dots, which is the central re¬ 
sult of this letter. The beamsplitter response contains a 
cross-polarisation contribution of 28±5 %, which is de¬ 
rived by integrating the measured beampatterns of the 
antenna 26 over the opening angle of the source (Supple¬ 
mentary Note 1). The other contribution to the error 
bars on the data is given by the uncertainty in the exact 
beamsplitter thickness 48 =b 2 /im. 

The power as a function of frequency that arrives at 
the detector waveguide input can be calculated using 

P(y) = 1 2 -J n A{u)B ^ Tbb) c(u)G(u, si)dn, (4) 
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FIG. 2. The measured phase response (dots) of the microwave 
resonator as a function of the frequency of the pair-breaking 
radiation and normalised to one. Error bars are shown every 
third point. The green line represents a calculation of the 
power absorption of the superconducting transmission line. 
The red line is a simulation of the pair-breaking efficiency 
(not normalised) that arises due to the different quasiparti¬ 
cle distributions at different excitation frequencies. The blue 
dashed line combines the two effects. The blue and green lines 
are both normalised to one. 


with c the speed of light, U the solid angle, A{y) the 
transmission of optical elements (filters, beamsplitter), 
B(v,Tbb) the brightness of the source given by Planck’s 
law, C(u) the antenna efficiency and G(z/, Q) the an¬ 
tenna gain pattern. The factor ( cjv ) 2 reflects a single 
mode throughput. For the purpose of the present exper¬ 
iment it is sufficient to know the relative power at each 
frequency. As discussed in Ref. 26 the beam patterns 
are measured in three frequency windows: 290-350 GHz, 
640-710 GHz and 790-910 GHz. The difference in the 
directivity for these bands is compensated by the differ¬ 
ence in the part of the source that they capture. The 
brightness of the blackbody at the measured frequen¬ 
cies can be well described in the Rayleigh-Jeans limit, 
where B(u,Tbb) = 2kTBB^ 2 /c 2 , which exactly compen¬ 
sates the frequency dependence due to the throughput 
(c/v ) 2 . The antenna efficiency is the only element from 
Eq. 4 that introduces a frequency dependence, as shown 
in Supplementary Figure 2. Using Eq. 4 we estimate 
the absorbed pair-breaking power from the FTS to be 
15 nW. The transmission of the optical elements A(v) 
is taken into account in the correction of the measured 
response as explained above. 

In Fig. 2, starting from 200 GHz, we observe no re¬ 
sponse until 320 GHz where the absorption rises dras¬ 
tically because photons have enough energy to break 
Cooper pairs (2A). This steep rise in response is partially 
the well-known absorption edge of the superconductor 28 : 


the frequency dependent absorption of a plain supercon¬ 
ducting film through the complex sheet impedance. How¬ 
ever in this experiment the antenna collects the radiation 
and launches it as a travelling wave into the MKID CPW. 
For frequencies well above the gap, it takes only 1 mm 
to absorb 90% of the power, thus ah power is absorbed 
in the detector volume. Therefore the non-monotonous 
sheet resistance for frequencies above the gap does not 
affect the measured response in this experiment, which 
is crucial to make the non-equilibrium response of the 
superconductor visible. The percentage of the power ab¬ 
sorbed in the CPW line is calculated using the atten¬ 
uation constant of a CPW 29,30 based on the frequency 
dependent complex conductivity of the Ta film following 
Mattis and Bardeen 8 . We assume the absorption length 
to be 10.4 mm, twice the length of the resonator. Ra¬ 
diation that is not absorbed (only for hv < 2A) will be 
reemitted by the antenna. The effective temperature of 
the superconductor is assumed to be 1.2 K, to be con¬ 
sistent with a linear response to 15 nW of FTS signal 
( N qp is dominated by the effective temperature). The re¬ 
sulting frequency dependent absorption is shown as the 
green line in Fig. 2. The maximum around 550 GHz is 
due to the simulated efficiency of the antenna, which is 
also taken into account (see Supplementary Fig. S4). 

For frequencies higher than 400 GHz, the power re¬ 
ceived by the antenna is fully absorbed in the detector 
waveguide. However, in Fig. 2a we observe a drop in 
the response close to 650 GHz (4A) by about a factor 
of two, after which the response increases again. Hav¬ 
ing taken into account ah frequency dependent power 
contributions, the only parameter left is the frequency 
dependence of the non-equilibrium response of the su¬ 
perconductor, represented by rjpb in Eq. 1. 

The non-equilibrium distribution of quasiparticles is 
calculated using a quasiparticle creation term that de¬ 
scribes the probability of creating a quasiparticle at 
a certain energy by breaking a Cooper pair following 
Eliashberg 5,31 . In steady state, the injection of quasipar¬ 
ticles at that energy is balanced by electron-phonon inter¬ 
action (scattering and recombination). The kinetic equa¬ 
tions for the non-equilibrium quasiparticle- and phonon 
energy distributions are solved following the approach by 
Chang and Scalapino 6 . The numerical procedure is ex¬ 
plained in Ref. 32. In the modelling we used an effective 
bath temperature of 1.2 K to account for the broadband 
power absorbed by the detector. It was not necessary for 
this temperature to be exact as rjpb is not strongly depen¬ 
dent on the bath temperature at low reduced tempera¬ 
tures (here T/T c = 0.27) 7 . The resulting distribution 
functions f(E) for constant absorbed power are shown 
for various frequencies in Fig. 3a. For higher excitation 
frequencies there are more quasiparticles with a higher 
energy, and therefore less weight close to the gap, where 
the resonator is sensitive. We therefore expect the max¬ 
imum resonator signal at v = 2A /h and a minimum at 
v — 4A /h. The pair breaking efficiency, i] p b, calculated 
from these distribution functions is shown in Fig. 2 (red 
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line). Combined with the frequency dependent absorp¬ 
tion (dashed blue line) it clearly describes the main shape 
of the measured response. 

It was shown by Guruswamy et al. 7 that the behaviour 
of /(E) for frequencies higher than v = 4A /h crucially 
depends on the phonon trapping factor. When phonons 
are released due to scattering or recombination, the ra¬ 
tio of their escape time r esc and the pair-breaking time 
Tpb determines how many quasiparticles can be generated 
from a single incoming photon. Only for T esc /r p b > 1 can 
Tjpb increase at energies above 4A. r p b is material depen¬ 
dent and equals 2.3xlO -11 s for Ta 33 (2.8xlO -10 s for 
Al). For 200 nm Ta on Si we obtain r esc = 2 ns 34 , which 
gives a trapping factor of 87, which makes Ta a favourable 
choice (over e.g. Al) to experimentally address the effect 
of phonon trapping. It is difficult to estimate the pre¬ 
cise trapping factor because the substrate is a relatively 
thin membrane and because of the Nb seed layer, but it 
is certainly large. In practice, r\ p b is the same for trap¬ 
ping factors of 15 and higher 7 . The minimum at 4A due 
to phonon trapping, which we observe in Fig. 2 quali¬ 
tatively distinguishes the non-equilibrium response from 
other frequency dependent phenomena. 

When Cooper pairs are broken, the created high en¬ 
ergy quasiparticles relax back to energies close to the gap 
on a timescale of 0.1-10 ns 20 . The response can well be 
described by an effective number of quasiparticles N qp 7 
using ijpb as in Eq. 1. From the calculated /(E), we 
derive rjpb , n qp (quasiparticle density) and 02 . The (al¬ 
most) linear relationship between those properties (Eq. 
1) is demonstrated in Figures. 3b and c. We emphasise 
that although these figures and Eq. 1 suggest a sim¬ 
ple effective N qp , the knowledge of the microscopic /(E) 
is needed to get that N qp at a certain (P, u) through 
Tjpb. An effective N qp or effective temperature model 35,36 
which may well describe a measurement at one specific 
frequency, even of the full /(E) 22 , would give constant 
N qp for constant absorbed power, regardless of the ex¬ 
citation frequency. It would therefore not explain our 
observations in Fig. 2. 

The qualitative agreement between measurement and 
simulation in Fig. 2 is very good, especially the peak 
around 2A and the characteristic 4A point are well rep¬ 
resented. The deviation that occurs at higher frequencies 
is most likely due to an incomplete understanding of the 
combination of the FTS system with the lens-antenna. 
Except for the mentioned uncertainties, the antenna effi¬ 
ciency and absorption length are not independently mea¬ 
sured and the removal of the residual ripple in the re¬ 
sponse is not exact. To get a deviation smaller than 
the 10-15% achieved now, one would need a complicated 
calibration with a bolometer with better sensitivity than 
the Ta MKID coupled to the same lens-antenna. We note 
that the characteristic impedance of the waveguide (Zq) 
is also frequency dependent, but it changes the power 
transmitted from antenna to waveguide by only 0.1%. 
The diffusion length of quasiparticles in Ta is 2 /1 m based 
on a recombination time of 50 ns (for an effective tem- 



FIG. 3. (a) The calculated quasiparticle distribution func¬ 
tions /(E) as a function of energy for different excitation en¬ 
ergies (2A, 3A, 4A, 5A and 6A). The lines for 5A and 6A 
coincide. The absorbed power was kept constant, and the re¬ 
sulting variation is thus only an effect of the frequency of the 
absorbed photons. The inset highlights the differences in /(E) 
on a small energy scale, (b) The change in the quasiparticle 
density (n qp ) and (c) the imaginary part of the complex con¬ 
ductivity 5(72/(J 2 as a function of the pair breaking efficiency 
r)pb • The different points are calculated at different excitation 
frequencies, 320 - 1100 GHz. The red lines are linear fits to 
the simulated points 


perature of 1.2 K) and a diffusion constant of 0.8 cm 2 s -1 
(Refs. 37 and 38 and the measured resistivity). The ef¬ 
fect of diffusion of quasiparticles from the central strip 
at the antenna feed is therefore negligible. Furthermore 
we checked that for this geometry radiation losses are a 
factor 10 lower than absorption in the superconductor 39 . 

The measured energy gap in the FTS response occurs 
at 324 GHz, corresponding to a T c of 4.4 K, assuming 
2 A = 3.52&#E C . This is consistent with the minimum 
response in Fig. 2 at around 650 GHz (4A). However the 
DC-measured T c of this film is 4.77 K, although most of 
our previous Ta films have also shown a T c of 4.4 K 40 . We 
presume that the Nb seed layer is thicker than anticipated 
giving a thin layer with a somewhat higher T c dominating 
the DC transport, whereas the radiation absorption is 
dominated by the lower gap in the thick Ta top layer. 

The microwave readout power can strongly affect the 
response of a microwave resonator 16,32 . In this experi¬ 
ment we can neglect effects due to the absorbed readout 
power (1.8 nW), which is much smaller than the absorbed 
pair-breaking signal (15 nW). Readout power effects are 
only expected in the opposite limit 41 , which is neverthe¬ 
less important to investigate in the future. The observed 
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agreement of the measurements with the model is encour¬ 
aging. At the same time it underlines the importance of 
understanding and controlling these parameters to opti¬ 
mise superconducting detectors. 
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SUPPLEMENTARY INFORMATION FOR: THE NON-EQUILIBRIUM RESPONSE OF A SUPERCONDUCTOR TO 
PAIR-BREAKING RADIATION MEASURED OVER A BROAD FREQUENCY BAND’ 

P.J. de Visser, S.J.C. Yates, T. Guruswamy, DJ. Goldie, S. Withington, A. Neto, N. Llombart, A.M. Baryshev, T.M. 
Klapwijk, and J.J.A. Baselmans 




FIG. SI. (left) The raw FTS response as a function of frequency together with the applied corrections for the transmission 
of the beamsplitter (both p- and s-polarisation contributions) and the low-pass filter. The error bars on the calculated beam 
splitter response are due to an uncertainty in the beam splitter thickness induced by the affixation of the beam splitter to a 
metal frame, (right) The corrected response from the left panel, shown here as the black line, contains a residual setup ripple 
due to the source, which we correct for as shown by the red line. The resulting response (blue line) is the data shown in Fig. 
2 of the main article. 


I. CALCULATION OF THE CROSS-POLARISATION CONTRIBUTION TO THE FTS RESPONSE 

The FTS response was measured without polarisers, which leaves the possibility of a cross-polar contribution in the 
response of the antenna. The beampatterns for co- and cross-polarisation were measured and reported in Ref. 26. As 
an example we show the measured beam patterns for 350 GHz in Fig. S3. We correct the cross-polar response for 
a small misalignment of the polariser by subtracting a 9% co-polar contribution, the result of which is shown in the 
right panel of Fig. S3. To calculate the cross/co-polar power ratio we first subtract the baseline level (-15 dB, due to 
noise) from all patterns. We subsequently integrate the beampatterns over an angle of ±4.5°, which corresponds to 
the opening angle of the broadband globar source as shown in Fig. S2. The resulting cross/co-polar response ratio 
is 0.27 for 350 GHz, 0.24 for 650 GHz and 0.33 for 850 GHz. To estimate this cross/co-polar ratio for the wideband 
FTS measurement we take the average of the result for the three bands and as the error bar the standard deviation 
of the three measurements: 0.28 =b 0.05. 
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FIG. S2. A photograph of the globar source taken with an infrared camera, demonstrating the 4.5 degree opening angle of the 
source. The solid lines in the centre are the 3 dB contours of the measured co-polarisation beampatterns at the frequencies 
indicated in the legend. 



FIG. S3, (left) Measured beampattern in co-polarisation in a narrow band around 350 GHz. (middle) Measured cross¬ 
polarisation pattern, (right) Measured cross-polarisation pattern with a contribution of 9% co-polarisation contribution sub¬ 
tracted. The color scales represent dB. 
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FIG. S4. (Left) The absorbed power in the superconducting coplanar waveguide detector as a function of frequency and as a 
fraction of the power at the CPW input. (Right) The efficiency of the antenna from a CST simulation, including impedance 
mismatch of the antenna to CPW, front-to-back ratio, spillover and lens reflection. The multiplication of the results in these 
two panels is the simulated green curve in Fig. 2 in the main text. 






